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(N| ■ ABSTRACT 

The Ka resonance complexes in oxygen ions O I - O VI are theoretically com- 
puted and resonance oscillator strengths and wavelengths are presented. The 
highly resolved photoionization cross sections, with rclativistic fine structure, 
CNJ ■ are computed in the coupled channel approximation using the Breit-Pauli 

R-matrix method. A number of strong Ka resonances are found to be appre- 

'. 

, ciable, with resonance oscillator strengths f r > 0.1. The Ka resonance wave- 

^ cj lengths of O I-O-VI lie in a relatively narrow wavelength range 22 — 23. 5A, and 

I ■ the X-ray opacity in this region should therefore be significantly affected by K 

O ■ 

' — > L transitions in oxygen. The results should be useful in the interpretation 

of soft X-ray spectra observed from Chandra and XMM-Newton . 
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1 INTRODUCTION 

K-shell excitation leads to the formation of resonance complexes in all atomic systems, with 
the exception of H- and He-like ions. The Ka complex in particular, associated with the Is 
— ► 2p transition, is expected to be generally strong. The same transitions in H- and He-like 
ions are well known and have been observed in many Chandra and XMM-Newton sources 
(e.g. Branduardi-Raymont etal. 2001, Lee etal. 2001, Ness etal. 2001, Sako etal. 2002); 
they are bound-bound transitions and their wavelengths and oscillator strengths have been 
determined to very high accuracy for most elements (e.g. Drake 1979, Savukov, Johnson, 
and Safronova 2002). However, there is very little reliable data in literature for resonances 
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since they involve autoionizing states, not readily computed theoretically or observed exper- 
imentally. 

We refer to lines as due to transitions between initial and final bound levels, and to reso- 
nances as due to transitions between a bound level and an interacting bound and continuum 
state (quasi-bound). The resonance may decay into the continuum with a finite lifetime, 
leading to a resonance or autoionization width in energy; it may also decay radiatively to a 
bound level as in the di-electronic recombination process. In astrophysical observations, with 
resolution usually lower than typical autoionization widths of the order of meV, line and 
resonance features look similar. In X-ray work, resonances are often referred to as 'inner- 
shell lines', although resonances need not be associated with inner-shells excitations alone. 
While the relevant transition matrix elements may be computed by atomic structure codes 
by explicitly coupling the bound and continuum wavefunctions in an isolated resonance ap- 
proximation, the well established close-coupling approximation from atomic collision theory 
enables a complete description of the coupled wavefunctions and interference effects in res- 
onances. Following traditional atomic physics usage, we refer to a complex as a set of lines 
or resonances belonging to the same principal quantum number. (Eissner and Nussbaumer 
1969, and references therein). 

Recently, theoretical calculations of wavelengths and resonance oscillator strengths for 
the O VI KLL resonances (Pradhan 2000) predicted a feature at 22.05 A that was later 
detected for the first time in the Chandra spectra of Seyfert 1 galaxy MCG— 6-30-15 by 
Lee etal. (2001). Similar calculations have now been carried out for the Li-like ions, C IV, 
O VI, and Fe XXIV (Nahar, Pradhan, and Zhang 2001). These resonance oscillator strengths 
are needed in order to interpret the X-ray spectra of many sources. Of prime interest is the 
plasma in active galactic nuclei, photoionised by and surrounding the central source. We may 
expect that all ionization states of oxygen might be present in such sources, and should be 
identifiable through their Ka resonance complexes, potentially leading to the determination 
of column densities and ionization fractions. Identifications of K- and L-shell transitions 
in Oxygen and several other elements have also been made in the XMM-Newton spectra 
of Mrk 766 (Ogle etal. 2002, Mason etal. 2002), MCG-6-30-15 6 (Branduardi-Raymont 
etal. 2001, Lee etal. 2001), and other sources (Kaastra etal. 2002, Blustin etal. 2002). To 
enable the identification and analysis of such features, in this brief Letter we present the Ka 
resonance strengths derived from large-scale relativistic close coupling calculations for the 
photoionization cross sections of all ions from O I to O VI. 
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2 THEORY AND COMPUTATIONS 

A procedure for computing resonance oscillator strengths has been described recently by 
Pradhan (2000). It is well known that we may relate the line oscillator strength and the 
bound-free photoionization cross section apj to the differential oscillator strength 



de 



2z lfline , 6 < I 

where I the ionization potential, z the ion charge, v the effective quantum number at 



e = — ^2 in Rydbergs, and a and ao are the fine structure constant and the Bohr radius 
respectively. The quantity j- describes the strength of photoabsorption per unit energy, in 
the discrete bound-bound region as well as the continuum bound-free region, continuously 
across the ionization threshold. 

The ^ reflects the same resonance structure as the opi in the bound-free continuum. 
Combining the two forms of ^ we therefore define, in the vicinity of a resonance, the inte- 
grated 'resonance oscillator strength' as: 

Wi — •//)=/ f d/(J '~* Jf) ) de = (—^) f anie; -h - J f )de, (2) 



I ae t \ de J \47r 2 aao / 
where Jj, J/ represent the initial bound and the continuum symmetries. Eq. (2) may be 
evaluated from the detailed api for the symmetries concerned provided the resonance profile 
is sufficiently well delineated. In practice this is often difficult and elaborate methods need 
to be employed to obtain accurate positions and profiles (the background and the peaks) 
of resonances. Furthermore, relativistic effects need to be included to differentiate the fine 
structure components. Using the coupled channel formulation based on the R-matrix and the 
relativistic Breit-Pauli R-matrix (BPRM) method (Burke etal. 1971, Berrington etal. 1995) 
a large number of photoionization cross sections have been calculated for all astrophysically 
abundant elements including resonance structures, particularly in the Opacity Project and 
the Iron Project works (Seaton etal. 1994, Hummer etal. 1993). We compute the ^ for all 
ions under consideration from photoionization cross sections from elaborate and extensive 
relativistic close coupling calculations using the BPRM method. In Eq. (2) the resonance 
oscillator strength f r is integrated over AE r , the energy range associated with the autoion- 
ization width(s). An advantage of the present method is that f r may be computed not only 
for an isolated resonance, but also for overlapping complex of resonances that may not be 
observationally resolved. 
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In the close-coupling or the coupled-channel approximation the bound or resonant (e + ion) 
system is represented by a wavefunction expansion over coupled levels of the 'core' or 'target' 
ion described by configuration-interaction type eigenfunctions. The presnet calculations for 
O-ions employ expansions for initial ground states and Ka resonances associated with all 
different (e + ion) continua (symmetry) using a number of configuration-interaction type 
target ion configurations, including those with Is- hole in the K-shell and 2s- hole in the 
L-shell. The extensive list of such spectroscopic configurations, and other computational 
details such as (e + ion) short-range correlation configurations, will be reported elsewhere. 
We note that these close-coupling calculations are based on the theory of atomic collisions, 
and are quite different from atomic structure calculations generally employed to compute 
line oscillator strengths. 

The resonant transitions from the ground state of an Oxygen ion into the Ka autoionizing 
resonances are: 
O I : ls 2 2s 2 2p 4 ( 3 P 2 ) + hv — ^ls2s 2 2p 5 ( 3 P ° 12 , 1 P°) 

O II : ls 2 2s 2 2p\*S° 3/2 ) + hv-^ls2s 2 2p^P 1/2<3/2/ , 5/2 , 2 D 3/2<5/2 2 S l/2 2 P 1/2>3/2 ) 

O III : ls 2 2s 2 2p 2 eP ) + hv^ls2s 2 2p 3 ( 5 S° 2 , s D° 2 3 , x D° 2 ? S{? P^, 1 P?) 

O IV : ls 2 2s 2 2p{ 2 P° /2 ) + hu-^ls2s 2 2p 2 ( 4 P 1/2t3/2/>5/2 2 D 3/2t5/2 2 S 1/2 2 P 1/2>3/2 ) 

O V : ls 2 2s 2 { l S ) + hv — >ls2s 2 2p{ 3 P° 12 , 1 P°) 

O VI : ls 2 2s{ 2 S l/2 ) + hv^\s2s2p{ 2 Pl l2 ^Pl lwl2 ). 

The resonances on the RHS of (3) autoionise into a free electron and a core ion in 
ground or excited states, consistent with total fine structure (e + ion) J7r-symmetries given 
as subscripts of the LS (( 2S ' +1 )L 7r ) terms. Whereas we locate the resonances due to all final 
continuum symmetries on the RHS, some are overlapping or are weak with little integrated 
resonance oscillator strength. 

The monochromatic opacity due to the Ka resonances may be obtained from the accu- 
rate bound-free photoionization cross sections that delineate autoionization shapes or pro- 
files. Fine structure also needs to be included explicitly since, as shown in this work, many 
fine structure components of resonances may be present. Because several coupled continuum 
channels contribute to each KLL resonance complex, as shown in Eq. (3), individual au- 
toionization widths are not computed (and may be overlapping); instead the whole complex 
is delineated as function of energy. Integration over the entire Ka photoabsorption range 
therefore yields the resonance oscillator strengths reported in the next section. 
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Calculations for low ionization states are more difficult than for higher ones since the 
electron-electron correlations are stronger relative to the dominant Coulomb potential, and 
many more electronic configurations, both for the ion eigenfunctions and for the (e + ion) 
system described by continuum and bound channels, need to be considered. 

Theoretical resonance spectroscopy, as formulated under the present approach, yields 
precise and detailed resonance positions, shapes, and oscillator strengths, compelmenting 
traditional line spectroscopy. 

3 RESULTS AND DISCUSSION 

Fig. 1 shows the Ka resonant complexes of all O-ions 1- 0V1. The dominant components 
and the peak positions of resonances in each ion are shown. As mentioned above, each 
complex has several components; their positions may or may not overlap depending on 
the exact energies of the contributing angular and spin symmetries. For example, O 1 and 
O II show one single peak because the contributing final Jn symmetries all lie at the same 
energy (n refers to the parity of the electronic state). On the other hand resonances begin 
to spread out with ion charge in O III and O IV, which also present more complicated 
resonance structures. The ground level of O III, ls 2 2s 2 2p 2 ( 3 Po), photoionizes into the J 
= (1)° continuum of ls2s 2 2p 3 with the three components shown. O V has only one final 
J7t( 1 P°) since the initial level is ls 2 2s 2 ( 1 5 , ). The KLL O VI is a twin component system 
discussed in detail in earlier works (Pradhan 2000, Nahar etal. 2001). However, the O VI 
resonance at 21.87 A is about an order of magnitude weaker than the stronger one at 22.05 
A ; the latter was identified in the MCG— 6-30-15 spectrum by Lee etal. (2001). Most of 
these results for O-ions have been obtained for the first time (prior works discussed below). 

We note that the tentative identification by Lee etal. (2001) of the O I resonance in 
MCG— 6-30-15 appears to be slightly higher in energy than the present value, and previous 
experimental and theoretical values at ~ 23.5 A discussed by Paerels etal. (2001), who used 
the non-relativistic theoretical cross sections for O I computed by McLaughlin and Kirby 
(1998) to analyse interstellar X-ray absorption (the theoretical model spectrum was shifted 
by 0.051 A to match the measured centroid wavelength of the O I resonance). Our O II 
Ka resonance is about 0.08 A lower than inferred from experimental data by Paerels etal. 
We also note that the computed wavelength for O V Ka resonance at 22.35 A is in good 
agreement with the rest wavelength of 22.334 A recently observed with XMM-Newton by 
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Blustin etal. (2002) from the Seyfert 1 galaxy NGC 3783. O III and O IV resonances in Fig. 
1 clearly illustrate the overlapping and asymmetric nature of autoionizing resonance profiles 
of the repective Ka complexes. 

In Table 1 we present the Ka resonance oscillator strengths for all O-ions, most with 
f r > 0.1. As evident from Fig. 1 several components are nearly degenerate. The f r given in 
Table 1 are therefore the sum over features lying approximately within O.OlA of each other. 
We also present autoionization 'equivalent widths' W a (meV), together with maximum peak 
values <J m ax, of the Ka resonances (note that f r = W a xo- max )- It needs to be emphasized 
that while the values in Table 1 may be utilised to treat Ka resonances on par with 'lines', a 
more precise treatment is to consider the entire photoionization cross section, up to energies 
encompassing the entire series of K-shell excitations, in calculating the bound-free X-ray 
opacity. Such calculations are more extensive and are in progress. 

Another potentially important consideration is the level population in excited fine struc- 
ture levels of the ground LS term. If the source plasma is moderately dense, N e > 10 6 cm~ 3 , 
then the excited levels may be significantly populated. Photoabsorption therefrom would 
lead to additional Ka features to augment those presented herein. These calculations are 
also in progress. 

4 CONCLUDING REMARKS 

The main conclusions of this work are as follows. 

1. The Ka resonance strengths of oxygen ions may be used to identify and analyse soft 
X-ray spectral features around 0.5 - 0.6 keV, given the Chandra and XMM-Newton resolution 
of ~ O.OlA. 

2. Although not spatially co-existent, O VI and O VII regions may both be present 
and detectable via K-shell X-ray transitions from a given source. A Li-like ion, such as 
O VI, is easily ionized, while the ionization of a He-like ion, such as O VII, is the most 
energetic of all ionization states of an element and may exist in a plasma over the widest 
temperature range (Pradhan 1982,1985). Thus the detection of O VI may be indicative of 
possibly a larger fraction of oxygen in O VII driven by recombination : (e + O VII) — > 
O VI. The O VI KLL doublet absorption resonances at AA 22.05 and 21.87 A (Pradhan 
2000) lie between the well known forbidden (T or 'z'), intercombination (T or 'x,y'), and 
allowed (V or 'w') Ka emission lines of O VII due to transitions 2( 3 S'i, 3 P 2 ° 1 , 1 P°) — > 
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l^So) at AA 22.101, 21.804, and 21.602 A respectively (e.g. Ness etal. 2001). To augment 
theoretical studies, new collisional calculations for electron impact excitation of O VII have 
been carried out (Delahaye and Pradhan 2002). In addition, self-consistent cross sections and 
rates O VI/O VII/O VIII have been computed (Nahar and Pradhan, in preparation) with 
the same eigenfunction expansion for both photoionization and recombination, and a unified 
(e + ion) recombination scheme including both radiative and dielectronic recombination in 
an ab initio manner using the BPRM method (e.g. Nahar and Pradhan 1992, Zhang etal. 
1999). The new O VII total and level-specific recombination rates and A- values, up to n = 
10 fine structure levels, should generally help the study of recombination lines, such as the 
O VII 'He/3, 7, 5, e' identified by Sako et al. (2002) in the XMM-Newton spectrum 

3. Analagous to the oxygen Ka resonances at ~ 0.55 keV (22-23 A ), there are complexes 
corresponding to other elements such as carbon at ~ 0.3 keV, nitrogen at ~ 0.4 keV, neon 
at ~ 0.9 keV. For example, at ~ 0.94 keV Lee etal. (2001) have identified a line from He-like 
Ne IX (Is 2 — ls2p). Calculations are in progress for these elements. 

Note to be added in proof: A recent experimental and theoretical study of Ka photoion- 
ization of O II (Kawatsura etal. 2002) reports the resonance oscillator strengths obtained 
from multi-configuration Dirac-Fock (MCDF) calculations. Their combined value for the 
transitions 4 S , 3 / 2 — 4 -P5/2, 4 -P3/2, 2 P1/2 is 0.765, compared to the present gf-value of 0.736 for 
the total Ka oscillator strength. It is expected that the other transitions not considered by 
Kawatsura et al. would be weaker, and if so then the present results should agree with theirs 
to about a few percent. They also report that the measured position is 1.5 eV lower than the 
MCDF value; it appears that the energies in Table 1 and Fig. 1 for O II may be higher by 
a similar amount, thus providing an estimate of the uncertainties in the resonance positions 
in the present work. 
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Table 1. Ka resonance oscillator strengths f r for oxygen ions 



Ion 


E r (Ryd) 


E r (KeV) 


MA) 


fr 


W (meV) 


a max (MB) 


O I 


38.8848 


0.5288 


23.45 


0.113 


31.88 


48.05 


O II 


39.1845 


0.5329 


23.27 


0.184 


23.21 


107.7 


O III 
O III 
O III 


39.5000 
39.6029 
39.7574 


0.5372 
0.5386 
0.5407 


23.08 
23.02 
22.93 


0.119 
0.102 
0.067 


27.00 
12.08 
24.27 


59.92 
114.5 
37.48 


O IV 
O IV 
O IV 


40.1324 
40.2184 
40.5991 


0.5458 
0.5470 
0.5521 


22.73 
22.67 
22.46 


0.132 
0.252 
0.027 


27.11 
14.32 
21.91 


66.15 
239.2 
17.00 


O V 


40.7826 


0.5546 


22.35 


0.565 


14.01 


549.0 


O VI 
O VI 


41.3456 
41.6912 


0.5623 
0.5670 


22.05 
21.87 


0.576 
0.061 


1.090 
12.16 


7142.0 
67.36 
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Figure 1. Photoionization cross sections of Ko resonance complexes in oxygen ions O I - O VI. 
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